Searches for very high energy gamma rays from blazars with CANGAROO-III telescope 

in 2005-2009 
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^ Abstract 

£3 .We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric 
& Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, 
performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in 
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the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra 
of these objects were obtained by analyzing Fermi/LAT archival data. Wide range (radio to VHE gamma-ray bands) spectral energy 
distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data, even though they are non- 
simultaneous, are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) 
t-h -radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 
are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of 
magnetic field and component size. 
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1. Introduction 

S3- 

j_j ■ Blazars are a sub-class of active galactic nuclei (AGN), with 
C3 spectral energy distributions (SEDs) that are characterized by 
double-peaked nonthermal emission which extends from radio 
to gamma rays. In addition, they are also characterized by rapid 
flux variability in all bands, polarization in radio and optical 
bands, and highly collimated relativistic jets. AGNs have been 
classified into various classes according to the observed fea- 
tures, which can be understood by differences of the angle be- 
tween the jet and our line of sight dominantly, and also circum- 
nuclear obscuration and relativistic beaming are candidates of 
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deciding factor [1]. Blazars include BL Lacertae objects (BL 
Lacs) and flat spectrum radio quasars (FSRQs), which have jets 
within a few degrees of the line of sight. Following the Whip- 
ple group's discovery of the first TeV gamma-ray blazar Mkn 
421 [2], blazars have became one of the most interesting class 
of object for VHE gamma-ray astronomy. 

Recently, the AGILE ^ and Fermi [01 gamma-ray space 
satellites have been launched, the Fermi/LAT detector has sen- 
sitivities about factor 30 enhanced from EGRET |5j|. Combin- 
ing simultaneous observational data from these satellites with 
VHE gamma-ray data from imaging atmospheric Cherenkov 
telescopes (IACTs) allow us to study multiwavelength behav- 
ior in a wide range. 

H 2356-309 is a high-frequency peaked BL Lac (HBL) at a 
redshift of z — 0.165 flg], which was discovered in X-rays by 
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the UHURU satellite Ym. Some authors classified this object 
as "extreme HBL" HI |9[] based on its synchrotron peak at ex- 
tremely high frequency. In 2004, the H.E.S.S. group discovered 
this object in the VHE gamma-ray band above 0.2 TeV at the 
10cr level of significance 1 111 , as predicted by Costamante 
et al. lfl2ll . The observed VHE gamma-ray spectrum between 
0.2 TeV and 1.3 TeV is consistent with a power-law with pho- 
ton index of T = 3.1. H.E.S.S. observations were additionally 
made from 2005 to 2007 and they found little flux variability 
on the time scale of a few years HI]. The CANGAROO-III 
telescope was pointed to H 2356-309 in 2005, but no evidence 
of VHE gamma-ray emissions above 750 GeV was found [ 14]. 
The flux upper limits are about one order of magnitude higher 
than the H.E.S.S. reported spectra. In the GeV gamma-ray 
band, Fermi detected this object with a low level flux 1 15, Hill , 
though EGRET and AGILE could not HUSH- 

PKS 2155-304 is a nearby HBL with a redshift of z = 0.1 16 
I2(ill . discovered in X-rays by the HEAO-1 satellite ELl. This 
object has been intensively studied in VHE gamma rays, since 
the Durham group first reported VHE gamma-ray emission 
lH. The CANGAROO group tried to find evidence of VHE 
gamma-ray emissions but without success (e.g., lEil l24l I25I0 . 
This object was confirmed as a VHE gamma-ray source by 
the H.E.S.S. group in 2005 Q. In July 2006, the H.E.S.S. 
group detected an extreme flaring of PKS 2155-304, and re- 
ported flux variation on the time scale of minutes [27], The flare 
triggered CANGAROO-III observations of PKS 2155-304 and 
VHE gamma rays above 660 GeV were detected at the 4.8cr 
significance level J3. During the same period, the MAGIC 
group also detected PKS 2155-304 above 400 GeV JH. Since 
then, many multi wavelength observations of PKS 2155-304 
have been carried out (e.g., i30l[3""lto . and part of the multiwave- 
length campaign in 2008 (August 25 to September 6) was re- 
ported in Aharonian et al. 132]. CANGAROO-III observations 
were also conducted during the 2008 summer, but no significant 
gamma-ray signal above 720 GeV was found lfl4Tl . The flux up- 
per limits were obtained in this observation slightly higher than 
H.E.S.S. observed spectra lf32h - This object was detected also 
in the GeV energy region by EGRET 
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PKS 0537-441, at a redshift of z = 0.894 |33|], is one of the 
most distant and the most luminous members of low-frequency 
peaked BL Lac (LBL) class, and may be part of a small galaxy 
cluster/group with several companions 113411 . Although clas- 
sic/normal BL Lacs have no features in their optical spectra, 
this source has several kinds of broad lines [35], and while the 
lines disappear in high states [36]. Therefore, PKS 0537-441 is 
discussed as a bridging object between the BL Lac and FSRQ 
classes (e.g., B37l 13810 . Although no detection of this source 
has been claimed in the VHE gamma-ray regime, GeV gamma- 
ray detections have been reported by EGRET iTTvL fTsh , AGILE 
fl^l . and Fermi lfl5l [l6ll at high significance levels, and GeV 
gamma-ray flares are often reported (e.g., J39L I40L 41, 4211. 



not [45, 46J. It was first detected in 2006 by the MAGIC group 
above 80 GeV [47], and became the most distant (identified) 
source of VHE gamma rays. The MAGIC group detection of 
the source again on 2007 January 16 [48], however, observa- 
tions by the H.E.S.S. group from January 18 to 21 did not de- 
tect this source [49]. At GeV energies, 3C 279 is a well known 
GeV gamma-ray blazar detected by EGRET USUI], AGILE 
fl^l . and Fermi 1 15, T^, and GeV gamma-ray flares have been 



reported several times (e.g., 11501,15 111 ). A multiwavelength cam 



3C 279 is a distant FSRQ, located at a redshift of z = 0.536 
1 43], and is one of the most intensively studied blazars over 
many years ll44Tl . In the VHE region, the Whipple group tried 
to detect VHE gamma rays several times in the 1980s, but could 



paign including Fermi/LAT was carried out from July 2008 to 
June 2009, and a change of the optical polarization associated 
with a GeV gamma-ray flare was reported ll52ll . During this pe- 
riod, the MAGIC group did not detect VHE gamma rays from 
this object 15311 . CANGAROO-III observations are also done 
during February and March 2009, and this results are reported 
here. 

In this paper, we report the results of VHE gamma-ray 
searches from these four blazars, which were observed with 
the CANGAROO-III telescope in the years from 2005 to 2009. 
This includes a reanalysis of H 2356-309 and PKS 2155-304 
data which were reported by Nishijima et al. lfl4ll . and new re- 
sults for PKS 0537-441 and 3C 279. We also discuss the mul- 
tiwavelength spectra using a simple one-zone leptonic model 
consisting of a synchrotron self-Compton (SSC) radiation and 
an external Compton (EC) radiation field. 

2. Instrument and observations 

The CANGAROO-III telescope system consists of four 
imaging atmospheric Cherenkov telescopes which are oper- 
ated near Woomera, South Australia (136°47'E, 31°06'S, 160m 
a.s.l.). These telescopes use altazimuth mounts and are placed 
at the east (Tl), west (T2), south (T3), and north (T4) corners 
of a diamond with sides of ~100 m. The Tl telescope, which 
is the oldest telescope, has not been operated for observation 
since 2004 due to its narrower field of view (FOV) and higher 
energy threshold. We also have not operated the T2 telescope 
since 2008 due to severe deterioration of its mirror. Only T3 
and T4 telescopes are used for analyses in this work. They 
have parabolic reflectors of 10 m diameter with an 8 m focal 
length, and each reflector is composed of 1 14 small spherical 
mirrors of 80 cm diameter and with an average radius of cur- 
vature of 16.4 m, which are made of Fiber Reinforced Plastic. 
Each imaging camera is an array of 427 PMT pixels of 0.17° 
size, with a total FOV of ~4°. The PMT signals were recorded 
by charge ADCs and multi-hit TDCs. Details of the mirror, the 
data acquisition (DAQ) system, the camera, the standard cali- 
bration method, and the total performance of the light collect- 
ing system are given in Kawachi et al. |54ll. Kubo et al. 115 511 . 
Kabuki et al. I5al . and Enomoto et al. 15711 . 

The CANGAROO-III observations of four blazars are sum- 
marized in Table Q] Each observation was made for 12-18 
moonless nights using the wobble mode, in which the pointing 
position of each telescope was shifted in declination by ±0.5° 
every 20 min from the target. The wobble mode enables si- 
multaneous observation of background regions with the target 
region. These observations, except for the case of 3C 279, were 
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basically carried out at zenith angles smaller than 30°. Details 
of each observation follow. 

The observations of H 2356-309 were carried out over 14 
nights during the period between 2005 July 8 and August 12 
with the three telescopes (T2, T3 and T4). Due to DAQ trouble 
with the T2 telescope, twofold coincidence data (i.e., T3 and 
T4) were used in the analysis. The total observation time is 
30.3 h and the mean zenith angle of the observations was 12.7°. 

The observations of PKS 2155-304 were made over 18 
nights from 2008 July 29 to September 27 with two telescopes 
(T3 and T4) for a total of 57.5 h, with a mean zenith angle for 
the observations of 14.6°. 

PKS 0537-441 was observed in 14 nights during the period 
between 2008 November 24 and December 29 with the T3 and 
T4 telescopes. The observations were carried out over 36.7 h 
and the mean zenith angle of the observations was 17.2°. 

The observations of 3C 279 were carried out over 12 nights 
during the period between 2009 February 24 and March 28 us- 
ing the twofold coincidence of T3-T4. The total observation 
time is 40. 1 h and the mean zenith angle of the observation was 
34.3°. Because the culmination altitude of the target position 
is low, the observations were carried out at zenith angles up to 
~ 60°. 



3. Data reduction and analysis 

3.1. CANGAROO-III 

Here, we briefly describe the analysis procedure for 
CANGAROO-III data, which is almost identical with that given 
by Kiuchi et al. 15811 . More details can be found in Enomoto et 
al. i57ll59h. After the CANGAROO-III standard calibration of 
the relative gain and the timing for each pixel [56], we have 
cleaned the images to eliminate night sky background photons 
and select shower images as follows. Each shower image is 
required to be a pixel cluster consisting at least five adjacent 
pixels, each of which has a signal larger than 5.0 p.e. with an ar- 
rival time within +30 ns of the average hit timing. Some of tele- 
scope performances depends on observational zenith angles, Z . 
In particular energy threshold of detected gamma rays becomes 
higher with the zenith angle rapidly when the zenith angle ex- 
ceeds ~30°. So we cut the data of 3C 279 with observation dur- 
ing 30° < Z . For other three objects, we use all data because 
more than 90% of the data were taken at the zenith angle less 
than 30°. In order to reject the low quality data due to the ef- 
fect of clouds and dew, any data taken during periods when the 
shower rate was lower than the typical rate for clear nights were 
not used. The peaks of shower rate distribution of each object 
are considered as the typical shower rate, which are 7.0, 6.8, 4.9 
and 6.7 Hz, for H 2356-309, PKS 2155-304, PKS 0537-441 
and 3C 279, respectively. We set shower rate thresholds to 
6.0 Hz for H 2356-309, PKS 2155-304 and 3C 279, and 4.5 Hz 
for PKS 0537-441. Furthermore, in order to ensure reliable 
arrival direction and energy estimations, i.e., to avoid the de- 
formation of the shower images by the camera edge, it was re- 
quired that none of the brightest 15 pixels of each image should 
be in the outermost layer of each camera. 



The moments of the shower images (length and width) are 
parameterized as defined by Hillas B6011 . The arrival direc- 
tions are reconstructed using those the intersection of image 
axes. After the event reconstruction, the FD value is calculated 
from width and length using Fisher Discriminant (FD) method 
1 61], which have a mathematically maximum separation be- 
tween the distributions of gamma-ray events and background 
hadron events. 

The number of detected gamma-ray events are estimated us- 
ing the FD fit method 15911 for each squared angular distance 
bin. As in our previous analysis H28T1 . to avoid the effect of 
position dependence of width and length in the FOV, each of 
background regions was employed in an area of the opposite di- 
rection relative to the center of the FOV. In this analysis, the FD 
distribution of background events is calculated for each ring- 
shaped region with area of 0.1 7r deg 2 , which radius is corre- 
sponded to the offset distance of gamma-ray analysis bin from 
the center of the target position. This method cancel out the 
position dependence of image parameters between the region 
of interest and the background region via the wobble motion, 
which switches each other region, and suppress the systematic 
difference. 

The energy threshold of the integral flux, E t h, is defined as 
the peak energy of Monte Carlo (MC) simulated gamma-ray 
events which survive pseudo hardware cuts and analysis cuts. 
For differential spectra, we use three energy bands correspond- 
ing to photoelectron bands of 0-80 p.e., 80-150 p.e., and above 
150 p.e., respectively. The mean energy of each photoelec- 
tron band is used as a representative energy of the band. The 
assumed VHE photon spectrum of each blazar for MC sim- 
ulation is a simple power-law with the following photon in- 
dices; F = 3.1 for H 2356-309, as determined by H.E.S.S. 
|Tll], r = 3.3 for PKS 2155-304, as reported by H.E.S.S. from 
their 2008 August to September observations [32], F = 6.0 
for PKS 0537-441, as a typical value for blazar at redshift of 
z ~ 0.9 with softening by extragalactic background light (EBL) 
photons, because there has been no detection in VHE region, 
and r = 4.1 for 3C 279, as reported by MAGIC JH. 

3.2. Fermi/LAT 

In order to obtain a wide-band high-energy SED, we use non- 
simultaneous GeV energy data taken with the LAT on the Fermi 
Gamma-ray Space Telescope [4]. We analyze the data between 
2008 August 5 and 2011 May 10, using the standard analysis 
package of ScienceTools-v9rl8p6. 

We selected the events with photon energies in the range of 
0.2-300 GeV. Two higher quality classes, the Diffuse class 
and the DataClean class, including both front and back events, 
were used, and the set of the instrument response functions 
of "P6_V3_DIFFUSE" was applied. We excluded events with 
zenith angles larger than 105° and time intervals when the rock- 
ing angle was larger than 52°, so as to reduce the contamination 
of Earth albedo gamma rays and the bright limb of the Earth. 

The standard unbinned likelihood analyses using gtlike were 
performed for our four regions of interest (ROIs), which are 
circular regions of 10° radius centered on the positions of each 
target blazars, H 2356-309, PKS 2155-304, PKS 0537-441, 
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Table 1 : Summary of CANGAROO-III observations for the four blazars 



Blazar 


Observation periods a 


Nights b 


^obs 




Target positions 6 




Year Begin End 




[his] 


[°] 


RA[° ] Dec[° ] 


H 2356-309 


2005 July 8 Aug. 12 


14 


30.3 


12.7 


359.783 -30.628 


PKS 2155-304 


2008 July 29 Sep. 27 


18 


57.5 


14.6 


329.717 -30.22 


PKS 0537-441 


2008 Nov. 24 Dec. 29 


14 


36.7 


17.2 


84.7098 -44.0858 


3C 279 


2009 Feb. 24 Mar. 28 


12 


40.1 


34.3(26.8/ 


194.046 -5.789 



(a): Details of observation periods, (b): The number of observation nights, (c): Each observation time, in units of hour, (d): The 
average observation zenith angle, in units of degrees, (e): Details of target positions in J2000 coordinate, (f): Averaged zenith angle 
of analyzed data after zenith cut of 6 Z < 30°. 



and 3C 279. We considered galactic and extragalactic diffuse 
emissions, and also considered the radiation from point sources 
which are cataloged as 1FGL [15] within the circle regions of 
15° radius from the ROI centers. The time variations of the 
sources within 10° from the center of the ROIs were taken into 
account, but for the sources outside that region the fluxes and 
the spectral indices are fixed at the values of the 1FGL catalog. 
In addition, we found several point-like sources near the tar- 
gets with a significance level of greater than at least 5cr, which 
are not listed in 1FGL catalog. We summarize these unknown 
source positions in Table|2] and these are considered as point 
sources in our analysis. Recently, second source catalog of 
Fermi/LAT (2FGL) is available [62]. Half of our additional 
sources might be appeared in the 2FGL catalog, but the oth- 
ers have no possible counterparts. This discrepancy might be 
caused by the different combinations of analyzed energy range 
and time range. 



4. Observational results 

After the data reduction described in previous section, the 9 2 
distributions of TeV gamma-ray candidate events are obtained 
by fitting the FD distribution of on-source data with those for 
MC gamma-rays and real background data. The left panels 
of Figs. HH3 show results for H 2356-309, PKS 2155-304, 
PKS 0537-441, and 3C 279, respectively. The number of ex- 
cess events in each O 2 bin are plotted with lcr statistical error 
bars. The hatched histogram represents our point spread func- 
tion normalized to the number of excess events in 9 2 < 0.06. 

The right panels of Figs. Q~H4]show the Fisher Discriminant 
(FD) distributions. Black filled circles with error bars are those 
for on-source events (0 2 < 0.06). Red solid and blue dashed 
histograms are the best-fit FD distributions for background 
events and for MC gamma-ray events, respectively. Green 
filled squares with error bars show the background-subtracted 
gamma-ray candidate events. 

In these figures, there is no significant gamma-ray excess 
from the direction of any of these blazars. We therefore cal- 
culated 2cr flux upper limits of gamma-ray emission for each 
object. The results are summarized in Table [3] In this table, 
the number of excess events, statistical significance, live time, 
assumed photon indices in MC simulation, estimated energy 
threshold, and the integral flux upper limits are listed. 



Meanwhile in the GeV gamma-ray bands, all of four targets 
are detected with high levels of significance. Our results of 
H 2356-309 show a very low flux level and a flat (r = 2.0) 
spectrum in the energy range of 0.9-70 GeV. The gamma-ray 
emission of PKS 2155-304 extends from 0.2 to 300 GeV with 
a hard (r = 1.9) slope. For PKS 0537-441 and 3C 279, 
curved gamma-ray spectra are detected, with gamma-ray emis- 
sion vanishing above energies of 140 GeV and 30 GeV, respec- 
tively. Our gamma-ray spectra are in good agreement with 
the 1FGL catalog HI. Time variations of PKS 2155-304, 
PKS 0537-441 and 3C 279 are seen in our analysis. Fermi/LAT 
GeV flux levels during CANGAROO-III observation periods 
show about 40% higher for PKS 2155-304 and 3C 279, and 
about 40% lower for PKS 0537-441 than long time averaged 
one. 

5. Spectral energy distributions and discussion 

Here, we discuss multiwavelength spectral energy distribu- 
tions (SEDs) including CANGAROO-III upper limits, GeV 
gamma-ray spectra of Fermi/LAT, and archival data, using a 
simple emission model. Our emission model is based on the 
one-zone homogeneous, time independent, leptonic jet model 
by Kataoka et al. 16311 . In this model, we take into account the 
synchrotron radiation of relativistic electrons for the low fre- 
quency component, and photons boosted by the same popula- 
tion of electrons via the inverse Compton scattering effect for 
the high frequency component. As the seed photons, both syn- 
chrotron photons and external thermal photons are considered. 
The former case is known as synchrotron self-Compton (SSC) 
radiation (e.g., I64L I65IH . and the latter as external Compton 
(EC) radiation (e.g., I66ll67ll). 

We assume there is a spherical blob, or jet component, of 
radius R with homogeneously tangled magnetic field B which 
propagates at relativistic speed with a bulk Lorentz factor r. 
The viewing angle 9, that is the angle between the jet orientation 
and our line of sight, is assumed to be the superluminal angle, 
in which case Doppler beaming factor 5 can be taken 5 - T. 
In the emitting blob, injection, radiative cooling, and escaping 
of electrons are assumed to be in equilibrium. We considered 
two models of a steady state electron energy spectrum. One is 
a simple power-law spectrum expressed as 
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Table 2: Summary of additional source position of Fermi/LAT analysis (RA [° ], Dec [° ]) 



H 2356-309 



PKS 2155-304 PKS 0537-441 



3C 279 



1 


(358.36, 


-30.58) a 


(327.97, 


-30.42) a 


(83.00, 


-48.46) a 


(192.74, -2.05) 


2 


(359.55, 


-28.89) 


(327.74, 


-27.69) a 


(83.01, 


-38.81) 3 


(198.14, -4.48) a 


3 


( 3.02, - 


-30.58) 


(325.52, 


-25.88) 


(79.67, 


-45.74) a 


(189.33,-4.39) 


4 






(335.70, 


-35.00) a 


(89.03, 


-43.92) a 


(187.93,-8.47) 


5 














(190.67,0.79) 



(a): These sources have a possible counterpart listed in 2FGL catalog. 
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Figure 1: B 2 distribution of excess events and Fisher Discriminant distribution for H 2356-309. (left): The plots are the number of excess events as a function of 
the squared angular distance d 2 . (right): Obtained Fisher Discriminant distributions within 6 1 < 0.06. The black filled circles are the observed events in the ON 
source region, the red solid and the blue dashed lines are the background and the gamma-ray components estimated by the FD fit procedure. The green filled squares 
indicate gamma-ray candidates, which are the excess events between observed events and estimated background events. 



Table 3: Summary of CANGAROO-III observation results 



Blazar 


Excess 3 


Significance b 


^live 


r d 


£th e 


/tj.l.(> £th)' 




[events] 




[his] 




[TeV] 


[cirT 2 s~'] 


H 2356-309 


-16.9 + 25.2 


-0.67 


23.3 


3.1 


0.68 


<2.6x 10 - 12 


PKS 2155-304 


59.5 + 29.9 


1.99 


41.8 


3.3 


0.60 


< 3.8 x 10~ 12 


PKS 0537-441 


20.8+18.0 


1.16 


26.6 


6.0 


0.51 


< 7.0 x 10~ 12 


3C 279 


-11.7 ± 12.9 


-0.90 


12.5 


4.1 


0.72 


< 1.7 x 10~ 12 



(a): The number of gamma-ray candidates and the errors, (b): Statistical significance levels of detected gamma-ray candidates, (c): 
Effective exposure times, (d): Photon indices of assumed power-law spectra in MC simulations, (e): Estimated energy thresholds 
of integral flux, (f): 2<x upper limits of integral flux above each energy threshold. 
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Figure 2: 6 1 distribution of excess events and Fisher Discriminant distribution for PKS 2155-304. The panels, including symbols and lines are same as Fig.^ The 
hatched histogram in the left panel represents our PSF normalized to the number of excess events. 
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and the other is a smoothed broken power-law expressed as fol- 
lows, 



JVe(y) = 



bi-k 



V 7brk I \ Tbrk / 



(2) 



where go is the normalization factor, y is the Lorentz factor of 
relativistic electrons, p\ and p2 are electron spectral indices in 
the range from y m j n = 1 .0 to y^k, and ybrk to y max , respectively. 
Due to the relativistic speed of the blob, the density of exter- 
nal photons injected from behind of the blob are greatly dimin- 
ished and the photons are redshifted in the blob frame. On the 
contrary, external photons injected from front of the blob are 
strongly boosted the density and the photons are blueshifted. 
Several kinds of seed photons of the EC radiation could be ex- 
pected, such as a nuclear continuum scattered by broad-line 
clouds, line emission from broad line clouds photoionized by 
the central source, and infrared photons from warm dust heated 
by the nucleus. As argued by Inoue & Takahara B85I1 . for simpli- 
fication of the EC radiation, we only suppose a single compo- 
nent of seed photons with a blackbody spectrum of temperature 
r ext , which is injected head on into the jet component. There- 
fore in our SSC model, the free parameters are B, R, 6), 
Go. Pi, Pi, Tbrk, and y max . If and when we consider the EC 
component as well as the SSC (here after SSC + EC model), 
two more free parameters, T ext , and L; so , are added, where Lj so 
is an injected photon luminosity in the rest frame of the object. 

Accurate estimation of the error of each physical parame- 
ter is very difficult in the case of large dimensional fitting. As 
one of indications, we select B, R and <5, which are thought to 
be closely related each others, and we estimate the 1 cr statisti- 
cal fitting errors in the three dimensional parameter space. The 
other parameters of the three dimensional space are fixed to the 



best derived ones. 

As is well known, VHE gamma rays are absorbed by Ex- 
tragalactic Background Light (EBL) photons due to e + e~ pair 
creation. There are various models of the EBL density (e.g., 
B68l l69L r70IO . In this paper, we use the EBL density model by 
Dommguez et al. 117111 . to estimate the amount of absorption of 
high energy gamma rays. 

5.1. H 2356-309 

Our SED for H 2356-309, including archival data, is shown 
in the Fig. [5] CANGAROO-III upper limits from observations 
in 2005 are indicated by the closed red circles, and the average 
Fermi/LAT spectrum obtained by the analysis of the data be- 
tween 2008 August and 2011 May is plotted by the open red 
squares. 



The average H.E.S.S. spectrum from 2004 to 2007 01311 . 
which includes our observation period, is indicated by the 
closed purple triangles. Our upper limits have not imposed a 
strict limit, however the GeV spectrum seems to connect to 
H.E.S.S. TeV flux smoothly. No time variation of either TeV 
gamma-ray flux lfl3ll or GeV gamma-ray flux from this object 
was found during these observational periods. 

We, then fit our SSC model to the observed spectra of gamma 
rays and X-rays, which include RXTE data observed simulta- 
neously to the H.E.S.S. 2004 observations [11]. The fit results 
are shown by solid lines in Fig. [5] where the lower component 
shows the synchrotron photons and the higher energy compo- 
nent comes from SSC photons. The observed data from X-rays 
to gamma rays are well reproduced. The best fitted parameters 
are listed in Table [4] 

The derived Doppler beaming factor, 6 = 59, is larger than 
the previously reported values (e.g., ITll[l3i[72ll ). and the mag- 
netic field strength of B - 0.012 G is much weaker than the 
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Figure 5: SEDs of H 2356-309. The red circles are 2<x upper limits from 
2005 July to August with CANGAROO-III. The open red squares are an av- 
erage Fermi/LAT spectrum. The upside-down purple triangles are the average 
H.E.S.S. spectrum of observations from 2004 to 2007 (lJ- The filled blue 
squares are the 2004 RXTE campaign spectrum, and the green triangles are 
archival data during the high state in June 1998, both are described in Aharo- 
nian et al. Hill and references therein. The light-gray crosses are archival data 
from the NASA/IPAC Extragalactic Database (NED). The bold magenta curve 
is the best fit to the observed SED by our SSC model. Two thin solid curves are 
the components of the synchrotron radiation and the SSC radiation. 



values so far reported (e.g., 73]). Fitting to the newly ob- 
tained very low flux and the fiat spectrum of the GeV gamma- 
ray band might causes these differences. However, derived pa- 
rameters must have large uncertainty due to loosely constrained 
synchrotron component. Furthermore, SSC model with a bro- 
ken power-law spectrum of electrons can be also reproduce the 
SED, but there is no significant improvement of the goodness 
of fit. So we employ the simplest model which has no break in 
the electron distribution. It may be able to determine whether 
the electron spectrum has a break or not if we have simultane- 
ous soft X-ray data. This model uncertainty is also one of large 
factors of parameter errors which we cannot take into account 
in here. 

Meanwhile, neither our SSC model nor SSC + EC model 
could reproduce the archival flux data of the radio and the opti- 
cal bands. The radio photons are considered to be emitted in a 
different region from the high energy emission, for instance, a 
inhomogeneous jet model 117211 is one such model. The optical 
emissions might be a combination of several components, such 
as a large amount of thermal emission from the host galaxy |@], 
the same component as the radio emission, and the synchrotron 
radiation from relativistic electrons. 

5.2. PKS 2155-304 

The wide-band SED of PKS 2155-304, observed at various 
times, is summarized in Fig. [6] The closed red circles indicate 
the upper limits calculated from the data taken in 2008 with 
CANGAROO-III. The Fermi/LAT data are also plotted by the 
gray open squares. Our flux upper limits in TeV energies are 
lower than the flux we observed in 2006 [28], and also lower 



than the reported fluxes by the H.E.S.S. team during multiwave- 
length campaign period in 2003 and 2008 which are plotted by 
the green triangles and purple inverted triangles, respectively. 
Although our observation period includes their 2008 campaign 
period, the results are not necessarily inconsistent with each 
other as PKS 2155-304 is one of the most variable blazars in 
the VHE gamma-ray region. 

Our radiation model is fitted to the data from the infrared to 
TeV gamma-ray energies, which includes archival infrared and 
optical data from NED and the simultaneous broad band spec- 
trum observed in the 2008 campaign period reported by Aha- 
ronian et al. i32ll . The report include Fermi/LAT results, and 
the GeV gamma-ray activity in this period is slightly higher 
than the average one which we obtained from the data between 
2008 and 201 1, though the TeV gamma-ray fluxes show a qui- 
escent state. Our SSC model with single power-law electron 
spectrum cannot reproduce the SED, and hard photon slope of 
soft X-ray band require the break of spectrum strongly. The 
solid line in Fig. |6]shows our SSC model prediction of the mul- 
tiwavelength SED assuming a broken power-law spectrum of 
electrons. The obtained physical parameters are summarized 
in Table |4] PKS 2155-304 has been well studied and a lot of 
effort was invested to understand its SEDparticularly by one- 
zone SSC model (e.g., III Ezl Q). Our model fits 
the synchrotron component very well, and the inverse Compton 
component is also well fitted. A large beaming factor 6 = 43 is 
obtained, which is a little bit larger than the previously reported 
values. Other parameters are in good agreement with already 
reported ones, although some of the models assume a double- 
broken power-law spectrum as an electron energy distribution. 

5.3. PKS 0537-441 

The observed SEDs of PKS 0537-441 at several epochs are 
compiled in Fig. [7] Our flux upper limits obtained from the 
observations in 2008 are the only data at TeV energies, which 
are plotted by the closed red circles in this figure, together with 
our analyzed GeV gamma-ray data of Fermi/LAT symbolized 
by the open gray squares. The quasi-simultaneous multiwave- 
length data observed in 2008 assembled by Abdo et al. 17911 . the 
data from AGILE and Swift, and archival data from NED are 
also plotted in the same figure. 

We apply our radiation model to fit the data taken in 2008, 
where our averaged GeV data were not used and the data from 
the campaign periods was used instead, because we have no 
data of radio to X-ray bands simultaneous to our CANGAROO- 
III observations. The simple one-zone SSC model fails to fit 
the observed SED well, so we take into account another com- 
ponent, the EC component, and the results are shown in Fig. 
|2l where the low-frequency component, mid-frequency com- 
ponent, and high-frequency component are the synchrotron ra- 
diation, SSC radiation, and EC radiation, respectively. From 
this figure, it is found that contributions from SSC photons and 
EC photons to the X-rays are roughly comparable, while the 
gamma-ray photons mainly originate from the EC mechanism. 
The physical values obtained by the fitting are listed in Table 
[4] Although the strength of the magnetic field obtained here is 
much weaker than the typical values reported before by others 
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Figure 6: SEDs of PKS 2155-304. The red circles are 2<x upper limits from 
2008 July to September with CANGAROO-III. The open squares are an aver- 
age Fermi/LAT spectrum. The purple inverted triangles are simultaneous spec- 
tra of H.E.S.S., Fermi, RXTE, and ATOM in 2008 The blue filled squares 
are the exceptional high state in 2006 July [31]. The green triangles indicate the 
2003 RXTE-H.E.S.S. campaign results [30]. The open diamonds are archival 
campaign data from 1994 November [78]. The light-gray crosses are archival 
data of NED. The orange crosses are also NED archival data but the data was 
used to SED fittings. The bold magenta curve is the best fit to the observed 
SED by our SSC model. Two thin solid curves indicate the synchrotron and the 
SSC radiation. 



(e.g., [35, 80, 81]), some models of this object proposed by Pu- 
cella et al. 18211 . which use a similar state of multiwavelength 
behavior with ours, predict a magnetic field as weak as the one 
we derived. From our model, the size of emitting blob is de- 
rived to be R = 4.1 x 10 16 cm, and the bulk Lorentz factor is 
obtained to be r — 24, which are consistent with other works 
1 35L [U, [82] . Although the best of the black body tempera- 
ture of EC seed photon is derived to be T ex[ = 2.3 x 10 3 K, it is 
difficult to infer the physical situations. Because this parameter 
is one of the most insensitive parameter in our model. 

It is seen from the figure that our SSC + EC model cannot 
reproduce the observed high radio fluxes which requires a huge 
size of the emission blob in our model and is implausible from 
the time scales of the observed flux variations 18311 . The bulk of 
the radio emission therefore probably originates much farther 
out along the jet. 

5.4. 3C279 

Many simultaneous multiwavelength campaigns of 3C 279 
have been performed to date. The SEDs obtained from some of 
those observations for wide ranges of frequencies are shown 
in Fig. [8] CANGAROO-III upper limits are plotted by the 
red closed circles. Open gray squares indicate GeV gamma- 
ray spectrum from our analysis of Fermi/LAT data. One 
of the simultaneous observation campaigns was conducted in 
February 2009, which data are symbolized by the purple in- 
verted triangles in this figure, just before our observations with 
CANGAROO-III. In this period, 3C 279 maintains an active 



Figure 7: SEDs of PKS 0537-441. The bold red circles with arrows are 2<x up- 
per limits in 2008 November and December with CANGAROO-III. The open 
squares are an average Fermi/LAT spectrum. The purple upside-down triangles 
correspond to the quasi-simultaneous radio to gamma-ray spectra in August to 
October of 2008 [72]. The blue closed squares are AGILE results in 2008 Oc- 
tober [82], and first AGILE catal og 11911 . The green triangles are Swift results 
during 2008 August to October [80]. The open diamonds are 2005 Swift and 
REM three simultaneous observations and archival data 18 ill . The light-gray 
crosses are archival data from NED. The best fitted curve to the SED by our 
SSC + EC model is indicated by the bold magenta curve. Individually, three 
thin solid curves are the synchrotron, the SSC, and the EC components, from 
the lower to the higher energy peaks, respectively. 



state in GeV band with flux levels clearly higher than the av- 
erage flux, while in the X-ray band it is at a relatively low flux 
level. 

For 3C 279, the SSC + EC model fits the observed SEDs 
from optical to gamma-ray bands well, where simultaneous 
data taken in campaign periods of 2009 are used for fitting. The 
model spectra are represented by the solid line in Fig. [8] and 
the best set of model parameters is summarized in Table [4] In 
this model, boosted EC photons are the main origin of gamma 
rays, and only a small fraction of SSC radiation contributes to 
the X-ray bands. 

Although in our model, the luminosity of the SSC compo- 
nent is relatively low compared to earlier works (e.g., [53, 80, 
84, 85il86[]), most of our parameters such as the size of moving 
blob R - 1.0 x 10 17 cm, the bulk Lorentz factor F = 19, and 
the strength of the magnetic field B = 0.56 G, are roughly con- 
sistent with them. The lowness of SSC luminosity is caused 
by a complex combination of the blob size R, the magnetic 
field B, and the electron density n e . The obtained seed pho- 
ton luminosity that is required from EC component photons is 
Lj so = 2.2 x 10 41 erg s _1 . If it is assumed for the accretion disk 
luminosity that the UV luminosity of L = 2 x 10 45 erg s _1 es- 
timated from the observations in the historical lowest state of 



3C 279 ISTJl, it could be probable that L iso = 2.2 x 10 41 erg 
s _1 corresponds to 0.01 % of the disk luminosity. The temper- 
ature of seed photons spectrum of EC radiation is derived to be 
Text - 5.7 x 10 3 K, but it is difficult to discuss its origin because 
this parameter is one of the most insensitive parameter to the 
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Figure 8: SEDs of 3C 279. The red circles are 2<x upper limits in 2009 February 
to March from CANGAROO-III. The open squares are an average spectrum of 
Fermi/LAT. The purple inverted triangles and green triangles are simultaneous 
multiwavelength spectra in the epochs of MJD 54880 to 54885 and MJD 54950 
to 54960, respectively [52]. The blue closed squares indicate a large gamma-ray 
flare state of 2006 February 23, and the open diamonds indicate the spectrum in 
the epoch of 2009 January 21 to February 1 l47ll53ll . The light-gray crosses are 
archival data from NED. The bold magenta curve is the best fit to the SED by 
our SSC + EC model. There are three thin solid curves, from the lower to the 
higher frequency components, the synchrotron, the SSC, and the EC radiation. 



SED. 

As in the case of the other three objects, the flux of the syn- 
chrotron component in this model is insufficient to explain the 
large observed radio flux. The possibility that the radio emis- 
sion regions in blazars might be much farther from the central 
engine has been discussed (e.g., [88, 8^]). 



6. Conclusions 

We have observed four selected blazars, H 2356-309, 
PKS 2155-304, PKS 0537-441 and 3C 279, with the 
CANGAROO-III imaging atmospheric Cherenkov telescope 
from 2005 to 2009. No statistically significant excess of events 
above 51 0-720 GeV from the direction of any of these objects 
was found, and we derived flux upper limits for VHE gamma- 
ray emissions. In addition, we analyzed GeV gamma-ray data 
between 0.2 GeV and 300 GeV taken with Fermi/LAT from Au- 
gust 2008 to May 2011. 

To derive some important physical parameters of these 
blazars, we consider a simple leptonic jet model to explain the 
multiwavelength SEDs including GeV and TeV spectra, even 
though non-simultaneous. The observed SED of H 2356-309 
(HBL) could be explained by a simple SSC model with a sin- 
gle power-law electron spectrum, and to keep a consistency 
with GeV spectrum, we need to assume a large beaming factor 
6 — 59 and weak magnetic field strength of 0.012 G. Radiation 
from PKS 2155-304, a nearby HBL, was well modeled by the 
SSC scenario, and obtained parameters are consistent with ear- 
lier works. PKS 0537-441, a luminous LBL, was studied and 



we found SSC + EC model could explain the observed multi- 
wavelength spectrum where the EC component is dominant in 
the gamma-ray photons. The SED of one of the distant FSRQ, 
3C 279, were also well explained by the SSC + EC model. 

Additionally, from our parameter fit results in Table |4] as the 
HBL to the FSRQ of blazar sub-classes, it is seen that the beam- 
ing factor becomes smaller, and in contrast the strength of the 
magnetic field becomes stronger and the size of blob becomes 
larger. Although we could not take into account of the uncer- 
tainties such as simultaneity of the data and the differences of 
models, these latter two tendencies are in agreement with the 
proposed blazar sequence (e.g., O90ll9lll '). 
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Table 4: Summary of the best parameter sets of the SED modelings 



[unit] H 2356-309 PKS 2155-304 PKS 0537-441 3C 279 

class HBL HBL LBL FSRQ 

~z 0165 0116 0894 053(5 

R [cm] (5.8!!^) x 10 15 (5.8+^) x 10 16 (4.1+^) xlO 16 (l.0!°;|) x 10 17 

B [G] (l.2+^) xlO- 2 (2.7+0; 7 ) x 10- 2 (4.7^) x 10" 1 (5.6^) x lO" 1 

r SQ+ 60 A ^+2-3 74+0.7 i n+0.8 

1 J -23 ^ -8.6 ^ -0.7 1 -0.8 
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[°] 
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1.3 


2.4 
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Go 


[cm _J y~ l ] 


1.5 x 10 4 


8.2 x 10-' 


1.5 x 10 4 


1.4 x 10 3 


pj 




-2.2 


-1.4 


-2.2 


-2.1 


P2 






-3.9 


-4.3 


-4.9 


Train 




1.0 


1.0 


1.0 
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Tbi-k 
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2.3 x 10 3 


1.2 x 10 3 


ymax 




8.1 x 10 5 


8.4 x 10 5 


8.7 x 10 4 


8.7 x 10 4 


^iso 


[ergs' 1 ] 
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2.2 x 10 41 


7"ext 


[K] 






2.3 x 10 3 


5.7 x 10 3 


Motel 




1.0 x 10 52 
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3.6 x 10 M 


5.3 x 10 s4 


«e 


[cm" 3 ] 
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2.0 x 10° 


1.2 x 10 4 


1.3 x 10 3 


e e 


[erg cm -3 ] 


5.7 x 10~ 2 


5.5 x 10- 4 


4.8x 10~ 2 


5.8 x 10~ 3 



The rows from upper to lower are as follows: blazar name, the class, z: the redshift, R: the size of radiative component with 1 cr 
statistical fitting error, B: the magnetic field strength with 1 cr statistical fitting error, F: the bulk Lorentz factor of the radiative 
component with 1 <x statistical fitting error, 6 V : the angle between the jet orientation and our line of sight, Qq: the normalization 
factor of electron spectrum, p\, pi', the electron spectral indices at low and high energies, y m j n , ■ybrk, y m!ix : the minimum, the break, 
and the maximum of electron Lorentz factors, L; so : the seed photon luminosity of EC radiation in blazar frame, r ext : the blackbody 
temperature of seed photon spectrum of EC radiation, N tota \: the total number of electrons, n e : the electron number density, and e e : 
the electron energy density, respectively. 
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